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We present a manifestly gauge-invariant description of Chern numbers associated with the 
Berry connection defined on a discretized Brillouin zone. It provides an efficient method of com- 
puting (spin) Hall conductances without specifying gauge-fixing conditions. We demonstrate 
that it correctly reproduces quantized Hall conductances even on a coarsely discretized Bril- 
louin zone. A gauge-dependent integer-valued field, which plays a key role in the formulation, 
is evaluated in several gauges. An extension to the non-Abelian Berry connection is also given. 
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Topological phase transitions have been of consid- 
erable interest in recent condensed matter physics. 1-3 
In lower dimensions, topological quantum numbers are 
known to play a crucial role in characterizing various 
phase transitions. A typical example is the integer quan- 
tum Hall transition, where quantized Hall conductances 
are given by Chern numbers associated with the Berry 
connection. 4-7 Its extension to the case of spin currents is 
also attracting much current interest. 8-11 These topolog- 
ical quantum numbers present a chance to characterize 
quantum liquids without using conventional symmetry 
breaking. 2,3 

Generically, the Chern numbers can be defined for 
quantum states with two periodic parameters. As shown 
below, they are given by an integral of fictitious mag- 
netic fields (field strengths of the Berry connection) over 
two-dimensional compact surfaces such as the Brillouin 
zone. In practical numerical calculations, however, we 
can diagonalize Hamiltonians only on a set of discrete 
points chosen appropriately within the surfaces. It is thus 
crucial to develop an efficient method of calculating the 
Chern numbers using restricted data of wave functions 
given only on such discrete points. In these calculations, 
a phase ambiguity of the wave function causes a gauge 
ambiguity for the Berry connection. Therefore, one must 
be careful if gauge-dependent quantities are used. 

In this letter, we propose an efficient method of cal- 
culating the Chern numbers on a discretized Brillouin 
zone. This is an application of a geometrical formulation 
of topological charges in lattice gauge theory. 12-16 We 
show that the Chern numbers thus obtained are mani- 
festly gauge-invariant and integer-valued even for a dis- 
cretized Brillouin zone. This implies that one can com- 
pute the Chern numbers using wave functions in any 
gauge or without specifying gauge fixing- conditions. For 
the purpose of demonstration, we apply our method to 
a simple model describing the integer Hall system. We 
find that even for coarsely discretized Brillouin zones, 
the method reproduces correct Chern numbers known so 
far. Our method can be useful in a practical computation 
for more complicated systems with a topological order 



for which a number of data points of the wave functions 
cannot easily be increased. 

To be specific, we focus on the Chern numbers in the 
quantum Hall effect. An extension to different topologi- 
cal ordered states is straightforward. The spin Hall con- 
ductances, for example, can be treated in a similar way. 
We consider a two-dimensional system in which the Bril- 
louin zone is defined by < fc M < 2n/q^ (/x = 1, 2 with 
some integers q^). Since the Hamiltonian H(k) is peri- 
odic in both directions, H{k\, k 2 ) = H(ki + 2n/qi, k 2 ) = 
H{k\, k 2 + 2tt /q 2 ), the (magnetic) Brillouin zone can be 
regarded as a two-dimensional torus T 2 . When the Fermi 
energy lies in a gap, the Hall conductance is given by 
o- xy = — (e 2 /h) J2 n °n, where c„ denotes the Chern num- 
ber of the nth Bloch band, and the sum over n is re- 
stricted to the bands below the Fermi energy. 4 ' 5 The 
Chern number assigned to the nth band is defined by 

Cn=^~. ! d 2 kF 12 {k), (1) 

ZTTl Jrp2 

where the Berry connection A fl (k) [fj, = 1, 2) and the 
associated field strength F 12 (k) are given by 4, 6:7 

A^k) = (n(fe)|9>(»), 

F 12 (k) = d x A 2 {k) - 8 2 A 1 (k), (2) 

with \n(k)) being a normalized wave function of the nth 
Bloch band such that H(k)\n(k)) = E n (k)\n(k)). In the 
above expressions, the derivative <9 M stands for d/dk^. 
We assume that there is no degeneracy for the nth 
state. 2 ' 3 The phase of the wave function is not yet de- 
termined here; that is, \n(k)) is defined on T 2 only up to 
its phase. 

If the gauge potential A^(k) is globally well defined 
over the continuum Brillouin zone T 2 , the Chern num- 
ber (1) vanishes because the torus has no boundary: It 
can be nonzero only when the gauge potential cannot 
be defined as a global function over T 2 . In this case, one 
covers T 2 by several coordinate patches and then, within 
each patch, one can take a gauge (that is, a phase conven- 
tion for the wave functions) such that the gauge poten- 
tial is a smooth and well defined function. In an overlap 
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between two patches, gauge potentials denned on each 
patch are related by a U(l) gauge transformation: 

\n(k)) -> e- 4A « \n(k)), A^k) - - i0 M A(fc). 

(3) 

The Chern number (1) is then given by a sum of the 
winding number of the U(l) gauge transformation along 
a boundary of a patch. As a consequence, the Chern num- 
ber is an integer. 

The above discussion is for the continuum Brillouin 
zone. Now suppose that we have data of wave functions 
only on discrete points within the Brillouin zone, as in 
actual numerical computations. A straightforward ap- 
proach for computing the Chern number (1) would be 
to replace all the derivatives by discrete differences and 
the integral by a summation. Namely one approximates 
the connection A fl (k)dk fl by 

A^k)dk^ = (n(fc)|5>(fc)>, (4) 

where 5^ is an infinitesimal difference operator defined by 
5 l _ L f(k) = f(k + 5kfj) — f(k) with 5k^ being an infinitesi- 
mal displacement vector in the direction fj, (its magnitude 
is \5kftl). Note that, to evaluate the difference, one must 
fix a local gauge with which the state |n(fc)) is smoothly 
differentiable near k. Under this local gauge, the field 
strength in the continuum is also approximated by 

F 12 (k)5k 1 8k 2 = [<M 2 (fc) - 5 2 A 1 {k)]5k 1 5k 2 . (5) 

Summing this Fi 2 {k)5k\5k 2 then gives the Chern number 
Cn in the limit \5k^ \ — > 0. However, this direct procedure 
can be costly in taking the limit if the Hamiltonian con- 
cerned is complicated. 

Here, we propose an alternative approach. Let us de- 
note lattice points kg (£ = 1, . . . , N\N 2 ) on the discrete 
Brillouin zone as 

k e = {k n , k j2 ), k^ = (j^ = 0, . . . , N„ - 1). 

(6) 

We assume that the state \n(ke)) is periodic on the lat- 
tice, \n(ki + N^fl)) = \n(ki)), where /t is a vector in the 
direction /i with the magnitude 2ir / (q^N^). Below, we 
set = q„N B (« ^ v) so that the unit plaquette is a 
square of the size 2w / (qiq 2 N B ) ■ 

We first define a U(l) link variable from the wave func- 
tions of the nth band as 

U^kt) = (n(k e )\n(k e + £)>/A/^), (7) 

where M^kg) = \(n(ki)\n(ke + /t)}|. The link variables 
are well defined as long as N^kg) ^ 0, which can always 
be assumed to be the case (one can avoid a singular- 
ity Af^(ke) — by an infinitesimal shift of the lattice). 
From the link variable (7), we next define a lattice field 
strength by 

F 12 (k e ) = lnC/ 1 (fc £ )C/ 2 (fc, + l)E/i(fo + 2)- 1 C/ 2 (fc £ )- 1 , 
-7T < -F 12 {k t ) < ir. (8) 

Note that the field strength is defined within the princi- 
pal branch of the logarithm specified in eq. (8). It is ob- 
vious that this field strength is invariant under the gauge 



transformation (3). Finally, we define the Chern number 
on the lattice which is associated to the nth band as 

5 " S ^-E^(fc£)- (9) 

First of all, we note that c n is manifestly gauge- 
invariant under eq. (3). This implies that we do not need 
to determine which gauge is adopted; any choice of gauge 
gives an identical number c„. Moreover, c„ is strictly an 
integer for arbitrary lattice spacings. To show this, we 
introduce a gauge potential 

A^k e ) - lnU„(k e ), -n < jA„(fc/) < w, (10) 

which is periodic on the lattice: A fl (ki + N^fi) — A^kf). 
Recalling definition (8), one finds 

F 12 (h) = AiA 2 (fc/) - AMh) + 2mn 12 (h), (11) 

where is the forward difference operator on the 
lattice, A^f(k e ) = f(k t + /t) - f(k e ), and n 12 (k e ) 
is an integer-valued field, which is chosen such that 
{l/i)Fi 2 (ki) takes a value within the principal branch. 
By definition, \n 12 (ki)\ < 2. From eqs. (9) and (11), we 
have 

c„ = ^ni 2 (fc £ ), (12) 

e 

which shows that the lattice Chern number c„ is an in- 
teger. 

The field strength on the lattice Fi 2 (ki) in eq. (8) re- 
duces to the one in the continuum Fi 2 {k)5k\8k 2 in the 
limit N B — > oo, where Sk^ = 2tt / (qiq 2 N B ) . Generically, 
the continuum field strength F 12 (k) has no singularity 
when the nth band is well separated from the neighbor- 
ing ones; that is, the energy gaps between them do not 
close, 

\E n (k) - E n±1 (k)\ ^0, (13) 

for any value of k G T 2 . This is the gap-opening condi- 
tion. 2 ' 3 One can expect, in general, that the problem is 
regular if the above gap-opening condition is satisfied. 
Then, the lattice field strength F\ 2 will be small enough 
for a sufficiently large N B and the lattice Chern number 
will approach the one in the continuum c„ — > c„ in the 
N B — > oo limit. Since both c„ and c„ are integers, we 
have c„ = c„ for N B > N B . The critical mesh size 
may be estimated by a breaking of the admissibility con- 
dition 12 - 16 

\F 12 (k e )\Ski5k 2 m |Fi2(fe/)| < for all k t . (14) 

It is expected that this N B is not very large for a stan- 
dard generic problem with the Chern number c„ « 0(1). 
Since the area of the Brillouin zone is 47r 2 / (qiq 2 ) , we can 
estimate the field strength as Fi 2 (kg) w ic n q\q 2 / (2n) . In 
this way, the critical mesh size is given by 

N B - « 0(^/2\c n \/( qi q 2 )). (15) 

That is, we can expect that our method reproduces cor- 
rect Chern numbers of the continuum even for a coarsely 
discretized Brillouin zone. This is another advantage of 
the present method. 
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As a function of U(l) link variables which satisfy the 
admissibility (14), the Chern number on the lattice c n is a 
constant function. To verify this, we note that a possible 
discontinuity of c„ as a function of link variables U^kg) 
occurs only when \Fi2^kg)\ = ir. Since c n is an integer 
which cannot continuously change, c„ remains the same 
as long as a configuration is smoothly varied under the 
admissibility (14). In other words, under the admissibil- 
ity, the space of U(l) link variables is divided into discon- 
nected sectors and the topological number c n is uniquely 
assigned to each sector. This is the basic idea behind the 
present construction. 12 The Chern number c„ is, more- 
over, a unique gauge-invariant topological integer which 
can be assigned to admissible U(l) link variables. 17, 18 

In the present context of the Berry connection, a 
gauge-invariant content of link variables is completely 
governed by the k dependence of the Hamiltonian H(k). 
Each of the topological ordered states with a nontriv- 
ial Chern number corresponds to the above nontrivial 
topological sector specified by the admissibility It is 
characterized by the lattice Chern number c n . In the 
continuum, on the other hand, the topological stabil- 
ity of the Chern number is assured by the gap-opening 
condition (13). The topological quantum phase transi- 
tions are thus characterized by the gap closing. Namely, 
nontrivial topological sectors of the continuum, each of 
which is a topological ordered state, are separated by 
the gaps. Correspondingly the Chern number of the to- 
tal bands, which is described by the non-Abelian Chern 
number, vanishes. 2,3 At the critical point at which the 
gap-opening condition breaks down, the field strength 
F\2{k) becomes singular at the gap-closing momentum. 
From a correspondence to the lattice case, we conclude 
that the admissibility condition cannot be satisfied by 
any finite JVb at that critical point. 

Our method can be extended to the case of the 
non-Abelian Berry connection A = ip'dif>, which is an 
M x M matrix-valued one-form associated with a mul- 
tiplet ip = (\ni), ■ • ■ , |^m))- 2 ' 3 ' 19 The associated Chern 
number is defined by cy, = J„ tr cL4/(27ri), an integral 
over a two-dimensional surface S with a generic (relaxed) 
gap-opening condition; E n (k) ^ E n '(k) for all k, where 
n E I and n! (£_ I for / = {n\, . . . ,njf}. 2,3 It turns out 
that the present lattice prescription is valid if one sub- 
stitutes the U(l) link variable by 

U^k e ) = -^—-det^(k e )^{k e + fi) (16) 

with the normalization constant N^k?) = 
| detip^ (kg)tp(kg + p,)\. We define the associated field 
strength and the Chern number on the lattice cy, by the 
same expressions as those for the Abelian case, eqs. (8) 
and (9). This C0 shares the features of the Chern number 
in the Abelian case c n . For regular problems, we have 
Cip = for a sufficiently fine discretization JVb > iVg. 

Having observed desired properties of our definition 
of the lattice Chern number, we now demonstrate how 
it works in a definite model. We consider the Hamilto- 
nian for spinless fermions in an external magnetic field: 
H = — t j) A e%di ^ c h where the nux P er plaquette 
on the coordinate lattice <p — 2~2n ^i-j / (^ 7r ) is p/l- 4 



For mutually prime integers p and q, the spectrum 
splits into q subbands. In the Landau gauge in the re- 
direction, the Hamiltonian in the fc-space is given by 
Hij(k) = -2t5ijCos(k y - 27r^j) - t(S i+ ij + 5 it j + i) - 
tS i+q -i, j e- iqk '' - t5 itj+ q-ie iqk * , (i, j = 1, q) with 
qi = q and q2 = 1. Bellow, we will present some results 
of applying our method to the middle subband of the 
= 1/3 (that is, q = 3) system. 




Fig. 1. (a) Field strength —iF\i(k) of middle band for <j> = 1/3 
system in 3 X 9 lattice Brillouin zone (JVb = 3), (b) the same 
result in 9 X 27 lattice (JV B = 9), and (c) field strength — iFi2(fc) 
of the same band in the continuum which is approximated by 
Im(5iA 2 - S 2 Ai) for \8ki\ = \5k 2 \ = 2vr/90. 

In Fig. 1(a), we show the lattice field strength F\2{k) 
in eq. (8). The (magnetic) Brillouin zone in the Landau 
gauge [0,27r/3) x [0, 2n) is discretized by 3 x 9 meshes. 
Note that the asymmetry of the Brillouin zone is simply 
due to the gauge choice; there is no x-y anisotropy in the 
present problem. The sum of F12 (k) over the mesh points 
gives c n — —2, which coincides with the known result for 
the present case. The same calculation but with 9 x 27 
meshes is shown in Fig. 1(b). It indeed gives the identi- 
cal Chern number c n = —2. Figure 1(c) shows the field 
strength F12 of the continuum. As expected, it is regular 
and of the order of unity. Comparing these figures, one 
can see that the field strength of the lattice system F\2 
converges to the one in the continuum Fyi up to a propor- 
tionality constant. Since the problem is regular, the field 
strength of the lattice system F\2 decreases as 2Vb in- 
creases. The admissibility is safely satisfied in Fig. 1(b) 
(JV B = 9), but JV B = 3 is close to the critical JVg (Note 
the scales in the figures.) 

It should be stressed again that the above lattice cal- 
culations can be performed in any gauge. We do not 
need specific gauge-fixing to make the gauge connection 
smooth. An arbitrary gauge (e.g., a phase choice of eigen- 
vectors given by a numerical library) can be adopted to 
compute the Chern number. 

As we have discussed, the lattice Chern number c n is 
closely related to the integer field n±2(ke) in eq. (11). 
To illustrate this point explicitly, we next plot the 
field ni2(kg). Since we must specify the gauge to do 
so (ni2(kg) itself is not gauge-invariant), we briefly de- 
scribe the method of gauge-fixing adopted here. 2,3 One 
first selects an arbitrary state \(f>) which is globally well 
defined over the whole Brillouin zone. Then the gauge 
can be specified by = P„|^)/JV* = \n) ■ {n\<j>) / N* , 
where P n = \n){n\ is a gauge-invariant projection and 
= |(<^|«)| is a gauge-invariant normalization which 
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ensures (n^\n^} — 1. A typical example of \(j>) is a con- 
stant state, but it can be a varying state as well. 

The integer fields ni2(ki) in several different gauges 
are depicted in Fig. 2, where the black and white cir- 
cles denote ni2 = — 1 and 1, respectively, whereas 
a blank implies ni2 = 0. It is clear that any of 
them gives the correct Chern number c„ = —2, 
that is, the number of black circles minus that of 
white ones is always two. The field ni 2 (k^) is gauge- 
dependent, but their sum is gauge- invariant. The fig- 
ures clearly show the gauge-independence of the lat- 
tice Chern number. In Figs. 2(a) and 2(b), we used 



(a) (2n/3,2n) (b) (2n/3,2n) ( c ) (2n/3,2«) M (2*3, 2k) 




(0,0) k (2n/3,0) (0,0) t (2rc/3,0) (00) i, (2tc/3,0) (0,0) ^ (2ji/3,0) 

"x K x *x K x 



Fig. 2. Configuration of integer field ni2(ki) in gauge speci- 
fied by state \<j>) (see text) over discretized Brillouin zones, (a) 
N B = 3, \<t>) = \<p Bl ) (b) N B = 3, \<j>) = |0 92 > (c) N B = 8 and 
\4>) = IHK/3.2)), and (d) N B = 8 and |0) = |^(J^ /4 . 2 )). Black 
(white) circles denote 7112 = —1 (1). 

the global gauges specified by the states \<fr gi ) — 
e i «( fc «'+ fc >')(l 1 -l,0) T and |0 S2 ) = e l< ?( fc - +fc «) (1, 1, 0) T , re- 
spectively. The term "global" means that the gauge- 
fixing condition |(</> ffi (A^)|n(&^))| ^ is satisfied at all 
lattice points kg. 

The meaning of the field ni2(ke) and the relationship 
between the present lattice formulation and the contin- 
uum one become much clearer by adopting a "patch- 
work gauge". To be specific, let us take a gauge con- 
vention specified by |</> 2 ) = (0, 1,0) T , but another gauge 
specified by |<fo) = (0,0, 1) T in some regions where 
the former is ill-defined. We show in Fig. 3 the am- 
plitude of the overlap between the wave function and 
trial states |</>2) and \<j>3). Figure 3(a) shows that the 
gauge specified by \(j) 2 } becomes ill-defined at two points 
near k 1 = (27r/3,7r/3) and k 2 = (tt/3, 4tt/3). There- 
fore, we first define, around these points, circular re- 
gions R r — {k | \k - fc 1 < r} U {k \k - k 2 \ < r} 
with an appropriate radius r, and we next check in Fig. 
3(b) that we can indeed take the second gauge specified 
by 1^3) safely in R r . This patchwork gauge choice is re- 
ferred to as the gauge specified by |0(i? r )). With this 
gauge, the wave functions \n^(ke)) and the correspond- 
ing gauge potential A^(ke) are smooth if lattice points 
are sufficiently fine. This implies that the integer field 
n\2 (kg) is vanishing within each region. Nonzero values 
of ni2(ke) are only allowed at plaquettes existing at the 
boundary of the regions. 



The field n\2{kt) computed with the above local 
gauge is shown in Figs. 2(c) and 2(d) for two differ- 
ent radiuses r. The figures clearly show that the inte- 
ger field ni2(ki) indeed acquires nonzero values only at 
boundaries of separated regions. The lattice field n\2(ki) 
carries information corresponding to the winding num- 
ber of the gauge transformation along the boundary of a 
patch in the continuum. 




Fig. 3. Amplitude of overlap [{<^j|n)| between wave function |n) 
and trial state \<j>i): (a) for trial state <j>2 = (0, 1, 0) T , and (b) for 
trial state $3 = (0, 0, 1) T . 

In this letter, we presented our method as an efficient 
technique for calculating the Chern numbers in an infi- 
nite system on the basis of a discretized Brillouin zone. In 
finite systems, however, the Brillouin zones are discrete 
from the onset. Therefore, the present method will also 
be useful for revealing topological orders of finite systems 
with possible many-body interactions. 2,3 
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